The Persistence Length of a Strongly Charged Rod-like Polyelectrolyte 

in Presence of Salt 
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The persistence length of a single, intrinsically rigid polyelectrolyte chain, above the Manning con- 
densation threshold is investigated theoretically in presence of added salt. Using a loop expansion 
method, the partition function is consistently calculated, taking into account corrections to mean- 
field theory. Within a mean-field approximation, the well-known results of Odijk, Skolnick and 
Fixman are reproduced. Beyond mean-field, it is found that density correlations between counte- 
rions and thermal fluctuations reduce the stiffness of the chain, indicating an effective attraction 
between monomers for highly charged chains and multivalent counterions. This attraction results 
in a possible mechanical instability (collapse), alluding to the phenomenon of DNA condensation. 
In addition, we find that more counterions condense on slightly bent conformations of the chain 
than predicted by the Manning model for the case of an infinite cylinder. Finally, our results are 
compared with previous models and experiments. 

PACS numbers: 61.25.Hq, 36.20.-r, 87.15.-v 



I. INTRODUCTION 

The behavior of charged polymers has received con- 
siderable attention since the early 70's. However, despite 
extensive efforts, much of the phenomena observed in sys- 
tems containing polyelectrolytes (PEs) is still not very 
well understood Q. PEs are frequently used in various 
industrial applications, such as stabilization of charged 
colloidal suspensions and in flocculation processes. They 
also are an essential ingredient of many biological sys- 
tems. These reasons motivated theoretical 0, @, 0, [| , 
experimental @, |§[ ^, |l^, and computer simulation 



jl2], [l3[ [IIJ [fj| investigations of PEs. For instance, DNA 
is known to be a particularly strongly charged polymer, 
bearing a charge density of one electron charge per 1.7A. 

Despite the strong electrostatic repulsion, PEs show 
a wide range of complex behavior, depending on the 
concentration of added salt and its valency. It is ob- 
served that with monovalent counterions, PEs are usually 
stretched and assume a rod-like conformation |l6|, [l7], [l8| . 
On the other hand, introduction of a small amount of 
multivalent counterions significantly reduces the rigidity 
of the chain |?], ||, [To| . Under certain conditions a PE 
may completely collapse into a globule-like conformation 
|il9[ |20|, For DNA, this is knows as DNA condensa- 
tion |23]. 

Even the problem of a single, uniformly charged, poly- 
mer in aqueous solution still poses a theoretical challenge 

Single- 
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chain models are much simpler than real experimental 
systems, as is the case for biopolymers in physiological 
conditions, or with synthetic PEs. However, much of the 



interesting phenomena characteristic to dilute PE solu- 
tions is still captured in the models despite the extended 
simplification. 

The first breakthrough in treating semi-flexible PEs 
analytically was made by Odijk f3(i|| , and independently 
by Skolnick and Fixman |5l| (OSF) by introducing the 
concept of an "electrostatic persistence length". The no- 
tion of persistence length 0, which measures correla- 
tions along the chain, is very useful in describing elastic 
properties of semi-flexible polymers in general, and PEs, 
in particular. According to the OSF theory, the total 
persistence length of the polymer can be written as a 
sum of two contributions: the bare persistence length, 
Iq, and an electrostatic one, l e . Electrostatic interactions 
are treated on a mean-field level and the charges on the 
polymer are assumed to be smeared uniformly. Within 
a linearized version of the Poisson-Boltzmann theory, 
the interaction between any two charges is screened and 
given (in units of k&T) by the Debye-Huckel expression, 
U DH (r) = Z B e" K 7r. The Bjerrum length l B = e 2 /ek B T 
is defined as the distance at which the electrostatic in- 
teraction between two ions of unit charge e equals the 
thermal energy k B T, where e is the dielectric constant of 
the medium. For water at room temperature, Ib — 7A. 
The inverse Debye-Huckel screening length is then k = 
[47tz(z+1)/bc] 1 / 2 . It depends on the concentration of salt 
c and the counterion valency z, where throughout this pa- 
per we will explicitly use a 1 : z salt. Finally, the polymer 
is considered to be intrinsically rigid, l ^> L. According 
to the OSF theory, the electrostatic persistence length 
then yields 



'p — h + 'osf — 
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where A = e/a denotes the average line charge density 
along the chain. 
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Since the OSF theory is strictly a mean-field theory, 
the effective interaction between charges o n th e polymer 
is always repulsive |32|, |33|. Indeed, Eq. (LI) indicates 
that the polymer becomes more rigid due to electrostatics 
because ?osf > 0. Experiments, however, clearly show 
that under some conditions, electrostatics may cause a 
reversed effect jj], p[ |9|, |22|, where enhanced chain flexi- 
bility results from a negative electrostatic contribution to 
the persistence length. In order to account for such be- 
havior, corrections to linearized Poisson-Boltzmann have 
been considered 
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|3J, |35|]. This was done in two 
steps. The first is to take into account the effect of coun- 
terion condensation. With strongly charged polymers, 
some of the counterions are loosely bound onto the chain 
and reduce the effective charge on the polymer. For a 
straight infinite cylinder, this is known as Manning con- 
densation |3|, |36| . The second step is to consider corre- 
lations between the ions and thermal fluctuations of the 
counterion density (3^, |3^] . Correlations between bound 
ions become more significant at lower temperatures. In 
the T — > limit, condensed ions are arranged on a pe- 
riodic lattice, similar to a Wigner crystal or a strongly 
correlated liquid . At higher temperatures such corre- 
lations are smeared out due to thermal fluctuations. The 
latter introduces another correction to mean-field theory 
coming from induced dipoles, similar to van der Waals 
interactions p5| . 

Both correlations and thermal fluctuations are mech- 
anisms that can cause the effective interaction between 
charges on the polymer to become attractive. Nguyen et 
al jM[ considered the former and calculated the persis- 
tence length of a polymer close to T — 0. Alternatively, 
Golestanian et al pa] considered the fluctuation mecha- 
nism. The two models do not agree with each other, and 
there is still no consensus on the question of which of the 
mechanisms is more significant at intermediate tempera- 
tures, which are used in experiments ^ [|| f^g 

The aim of the present work is to propose a model 
which takes into account both correlations and ther- 
mal fluctuations of a single, intrinsically rigid, charged 
polymer, immersed in a bulk and continuous dielectric 
medium. This will facilitate a closer and consistent exam- 
ination of the different mechanisms that cause the funda- 
mentally repulsive electrostatic repulsion to become at- 
tractive. Since the deviation from mean-field predictions, 
as seen both in experiments and in previous theoretical 
works, is so pronounced, an analytical understanding of 
the problem will be of high value. 

In the next section, we introduce our model for treat- 
ing a single, rod-like PE in presence of added electrolyte 
(salt) [E2|. In Section III, a mean-field approximation is 



used, and the well known result of OSF is reproduced. 
Section IV finds the first-order correction to mean field, 
taking into account both correlations and thermal fluc- 
tuations. This is the main result of the paper where 
a new expression for the electrostatic persistence length 
is obtained. This expression accounts for the observed 
attraction between monomers for strongly charged poly- 



mers and multivalent counterions. In the last sections, 
our results are compared to experiments and other theo- 
retical models. Our findings are further discussed in view 
of these comparisons. 



II. THE THEORETICAL MODEL 

Consider a polymer chain consisting of iV > 1 
monomers of length a each. Taking a worm-like chain 
approach the polymer is modeled as a spatial curve 
R(s), < s < L — Na, with a total persistence length 
Charges on the chain are assumed to be smeared with a 
positive constant line charge density X — e/a (one unit 
charge e per monomer size o) , while mobile ions are taken 
to be point-like charges. The system is immersed in a 
continuous dielectric medium with a dielectric constant e. 
For simplicity, we have assumed that co-ions are monova- 
lent, while counterions are multivalent, carrying a charge 
—ze. Namely, the chain is embedded in a 1 : z electrolyte 
solution. In order to account for the effect of counterion 
(Manning) condensation we follow the two-phase model 
introduced by Oosawa ||. The first phase is a ID gas of 
counterions that are bound to the polymer and can move 
only along the chain. The positions of / bound counteri- 
ons are denoted as R(si) . . . R(s/). The second phase is 
composed of free counterions in solution, in equilibrium 
with the ID gas. Since the Manning-Oosawa model re- 
gards the PE chains as an infinite-long cylinder, we will 
restrict ourselves hereafter to rod-like polymers which 
satisfy, ^ p ^ L. Finally, we will assume that the effect of 
the free ions is to screen out all electrostatic interactions 
[fTs] , [43| so that the interaction between any two charges 
(smeared charges on the polymer and bound z-valent 
counterions) is given by the screened Debye-Hiickel inter- 
action: F DH (r) = ZiZjl-Be~ Kr /r, where m and Zj are the 
valencies of the two respective ions. Because we employ 
a continuum approach, it is necessary to have < 1 (a 
is comparable to the size of the smallest molecule in the 
system - free monovalent ions and monomers) . Further- 
more, we require that the chain is long enough so that 
its contour length, L, is much longer than the screening 
length, kL ^> 1. These limits, a <C -C L, usually 
hold in experimental and physiological conditions, and 
are necessary conditions for our model. In particular, 
the no-added salt limit (k — > 0) is not covered by our 
model. 

We can proceed by writing down the grand-canonical 
partition function of the system. Up to a normalization 
constant it is 



Z = 



( IL1 ) 

where the path integral is a sum over all possible spa- 
tial conformations of the chain, /j, is the chemical poten- 
tial of the ID gas of bound counterions and is related to 
the counterion concentration in the bulk reservoir, Hq is 
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the Hamiltonian of a neutral chain with bare persistence 
length Iq, and H ln t is the electrostatic interaction Hamil- 
tonian. It consists of screened electrostatic interactions 
between all charged monomers and bound counterions 
and is written as a sum of three different contributions: 



H bh + H mh , where 

L L 



1 1 

2^ 



dsds' F UH (R(s)-R(s')), 







nib 



= JeE^ dh ( r (^)- r (^))- ( IL2 ) 



i=i j=\ 



All energies are dimensionless and written in terms of 
the the rmal energy fc^T. In this form, the integrations 
of Eq. ([1. 2) diverge as the terms contain also self inter- 
actions (for instance, when s — > s' in -ff mm ). All integra- 
tions, therefore, should have a lower cut-off at a distance 
of order a. 

In order to treat the interaction term analytically, it is 
more convenient to use continuous volume concentrations 
defined in the following way Q E3, [ifj, |47| 

m (r) = - / d.s(5(r-R(s)) 
a Jo 
I 

b (r) = £a(r- R(*)), (H.3) 

»=i 

where <j) m and b are the monomer and bound counte- 
rion concentrations at location r, respectively. These can 
be substituted into the partition function, Eq. ( II. 1 ) , by 
making use of the identity operator that couples discrete 
and continuous concentrations. This is done using the 
path integral representation of the Dirac delta function 



V<j> m (r) S m (r) - i J ds5{r - R(s)) 
o 

0^ m (r)DV m (r)exp|-i J d 3 r V> m (r) 
m (r) - i J ds5{r - R(s)) I (II.4) 



1 = J V<j) h (r)5 



(r)-£*(r-R(*,)) 

i=l 

y OV b (r)^ b (r)exp|-i ^ d 3 r ^ b (r) 
i 



The extra complexity of this method is the introduc- 
tion of two new auxiliary fields, denoted ip m and ip h , 
which couple to <fi m and 4> h , respectively. Substituting 
Eqs. (II. 4) and (11.3) into Eq. (II. 1), the partition func- 



tion reads 



Z = J VR(s) J] 2?0W&[R] 

V i— m,b / 



x exp(— H c 



6n[R] = exp[-iT id +- / d s V m (R(s))" 
a Jo 



+ i / d 3 r<?(> m (r)V> m (r) 



&[R] = ex Py d 3 r\an e^ (r) f(r) 
+l b (r)^ b (r)]} 



-f^cont 



d 3 rdV *(r)Z$(r')^ DH (r - r') 



— z z 



where we have defined 

n = eP/L. 
and £b was simplified in the following way: 



.j=o ' i=l 



(II.5) 



(II.6) 




x exp 



i J d 3 r^ b (r) |E (5 ( r - R (^)) 



.7=0 ' i=l J n / 



i£V> b (R(si)) 



exp 



exp 



n Q / d S e«/' b ( R ( s » 



an / d 3 re^ b(r V m (r) 



(II.7) 



It is easily seen that carrying out the integrations over the 
new fields tfi m , ip m , (j) h , ip h , reproduces Eq. (II. 1) exactly. 
However, the form of the continuous partition function, 
Eq. ( ]II.5| ), is better organized: single-body interactions of 
the monomer concentration <^ m and the bound counterion 
concentration </> b are collected into the terms £ m and £b- 
The two-body interaction term H coat has the form of a 
quadratic interaction between the concentrations vector 
field <fr, where the interaction between the two vector 
fields 3>(r) and <fr(r') is given by the matrix ZV DH (r—r'). 
In the above equation, we use vector notation to write the 
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interaction between the different fields as: 

*(r)Z#(r')F DH (r-r') 
= [0 m (r)0 m (r')-2z0 m (r)0 b (r') 
+ z 2 ct) h {r)(t) h (r')] F DH (r-r'). 



(II. 



This method can be easily generalized in order to account 
for any additional species the system may contain, or 
to different types of (non-electrostatic) interactions. For 
instance, a local interaction can be added to the matrix 
elements of ZV BR (r — r') in order to include excluded 
volume interactions. If counterions are replaced by more 
complex charged amphiphilcs, a hydrophobic attraction 
between the species can be added in a similar manner: 



1 -z 



V DH {r-r') + 



*anm ^mb 
V m b Vbb 



S(r- 



where (i,j=m or b) denote second viral coefficients 



Thus far, the partition function Eq. ( II. 5 ) is exact up 
to the general assumptions of the model — worm-like 
polymer, smeared charges on the chain, separation into 
two phases and screening by free ions. However, the in- 
tegrations cannot be carried out analytically and some 
approximations have to be made. The first-order correc- 
tion just reproduces the well-known mean-field results, 
as will be shown below in Sec. |lll| . Higher-order terms 
in the expansion represent corrections to mean field and 
will be presented in Sec. IV. 



where Z± is the partition function up to first order in 
if - The Fourier transform of </> b (r), ^> b (r) an d </> m (r), is 
denoted by i^V'k anc ^ ^k' respectively. The interaction 
Hamiltonian Hi consists of three contributions: the first 
term is the two-body Debye-Hiickel interactions. The 
second term is a bi-linear coupling of the concentrations 
fields, 4> m and </> b , with the auxiliary field ip h , generated 
by the bound counterions. As for the third term, we will 
later show that it is only a constant. The integrations 
over the degrees of freedom of the bound ions {<fi h 7 ip h } 
can now be carried out. 



Zi 



H C S,1 



PR(s)P0 m PV m Cme 



d 3 k 

(2tt) 3 



-(l- a n z) 2 V^ k DH CC k 



-an C(27r) 3 5(k)] 



(III.3) 



The effective interaction Hamiltonian between the 
monomers, i? c ff,i, is correct up to first order in ip h . 
The first term of H c g^ consists of screened electrostatic 
interaction (I^ DH ) between Fourier components of the 
monomer concentration field </> m . The interaction in- 
cludes a reduced charge density of (1 — an^z). Note that 
H c s,i depends on the conformatio n of the polymer R(s) 
through the definition of 4> m , Eq. ( II. 3 ). 



A. Averages and Correlations 



III. THE MEAN-FIELD APPROXIMATION 

The approximation method we use is a systematic expan- 
sion in powers of the auxiliary field ip h , similar to loop 
expansion in field theory Q, Q . Expanding to first or- 
der in ip h results in a mean-field approximation. The 
partition function takes into account the average interac- 
tion between the monomers and the bound counterions 
p4j . Note that this method is somewhat different than 
calculating the zeroth-order saddle-point approximation 
of the integral over UEL(s) in Eq. (11.5). For a detailed 



comparison between the two methods see Ref. J45|] . 

Expanding to first order in ip h , appearing only in £b 
we obtain 



~ expjy d 3 r[a7i o (l + ^ b (r))0 m (r) 

+*0 b (r)V b (r)]} (III.l) 
Applying a Fourier transform the partition function reads 



Zi = 



Hi = 



\i—m.b J 



d 3 k 



§*3 c i*-k^ k DH + i^ 1 



(2tt) 3 

ian Q ^\ - an ^(2TT) 3 6(k)] , (III.2) 



After presenting the partition function, our aim is to 
integrate out the degrees of freedom of bound ions and 
obtain averaged quantities up to first order in ijj h . Aver- 
ages over configurations of the ID gas of bound ions are 
defined as 



Z\ — Tr{ R 0m ^m}£ m Tr^b ^bje 



-Hi 



(Oh = 



T r {H.,<£»,V> m }£m e 



-He 



Tr, 



-Hi 



Tr { ^b } Oe- gl 

o— H e ft.-\ 



{.III.4) 



Taking derivatives of the partition function, different av- 
erages and correlation functions can be calculated. For 
instance: 

(^ b > 1 = iz(l - an z)V» H ^ 
<«>! = «>!«>! -(27r)V^ H ( 5(k 1 +k 2 ) 

<«>X = 

<«>i = <^ 1 ) 1 «) 1 -^(27r) 3 <5(k 1 +k 2 ) (III.5) 
It is interesting to observe that the effective interaction 



Hamiltonian -ff c ff.i of Eq. ([11.3) can be rewritten as 
d 3 k 



H 



eff.l 



±(l-an z)(^) 1 ^ 



(2tt) 3 
an <j)^{2TrfS(k)] . 



(111.6) 
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Up to a constant prefactor, the first term of H e ff,i is 
the interaction of the monomer concentration field m 
with the auxiliary field ip h averaged within mean-field 
approximation. Therefore, it is the averaged interaction 
between each monomer and the auxiliary field produced 
by the countcrions. 



an o [ f d3k >rl 



d J r / - - e ZK r - / dse lk R(s) 
L J r J (2^)3 6 a J 





^ /d 3 k / dse jk R ^5(k) 



B. Density of bound ions 

In the Manning - Oosawa's model j}[ [5(|, the PE is 
considered as an infinite charged cylinder, and the aver- 
age number of bound counterions of valency z per unit 
length is 



1 

zl B 



1 



z 2 h 



(III.7) 



where the dimensionless parameter q — zl^/a is used in 
Eq. 



III.7 



Condensation of z-valent counterions occurs 
for q > 1 and effectively lowers the value of q to unity, 
QeS = 1 Hi - We are more interested in the case when 
some of the counterions are condensed on the polymer. 
Note that bellow the Manning condensation threshold 
(q < 1) one should simply set n M = 0. 

The average number of bound countcrions per unit 
length is 



fc [R] = ~ / d 3 r0 b (r). (III. 



where the dependence on the polymer c onfor mation, 
R(s), is through the definition of b , Eq. (1L3). Since 
we treat here only single polymer chains, the system is 
assumed to be infinitely dilute, in the sense that each PE 
chain occupies only a small fraction of the overall system 
volume. Small changes in the polymer conformation are 
not expected to change the chemical potential of the free 
counterion gas which occupies the entire volume. As the 
ID phase of bound counterions is in equilibrium with the 
free counterion phase, their chemical potential is equal. 
We conclude that /x, and consequently n — e^/L, should 
not depend on the conformation of the chain in the dilute 
polymer limit. Of particular interest is the straight-rod 
conformation. In this conformation, the density of bound 
counterions should be consistent with the Manning the- 
ory, and n tot [R] should, therefore, satisfy 



'[rod] = 



,M 



(III.9) 



where n is the Manning value given in Eq. (III. 7) 



Within the mean-field approximation, the average num- 
ber of bound counter ions per unit l ength can be obt ained 
by substituting Eq. (flf) and Eq. flll.SD into Eq. ([III.8D 



[Rl 



d 3 r(^ b (r)) 1 



d 3 k 

(2tt) 3 



e ano0 k 



(111.10) 



1 



no / ds = uq. 



Equation (III. 10) gives us the connection between the 



average density of the bound counterions, n [R], and 
the chemical potential, \i (th rough no), up to first order 



in ip . According to (111.10), n^ ot [R] does not depend 



on the conformation of the polymer and is just equal to 
the concentration no. In the following section we will 
see t hat th is is s trictly a mean-field result. Comparing 



Eqs. fllH.9D and ( |III.10| ) we find that 



n 



C. The Persistence Length 

Substituting n a into the effective interaction, 
Eq. (pO), yields 



DH jja j,m 



L 



2z 2 



- X [R]-Ln M , 



(III.ll) 



where 



x[R] 



ds 





-«|R(s)-R(s')l 
ds'L-^ (III. 12) 



|R(,s)-R( s ')| 



Up to a constant, the integrand is just the screened 
Coulomb interaction between any two monomers. The 
integral is a sum over all such monomer pairs along the 
chain with charge density set at the Manning value (q = 1 
is equivalent to A — e/a = e/zls)- This is exactly the 
interaction Hamiltonian assumed by OSF for a polymer 
carrying a uniform line charge density of A = ej zIb- 

Odijk's method for calculating the electrostatic per- 
sistence length |30|] assumes small, constant curvature 
deformations from the straight rod conformation. The 
persistence length is then obtained from the rigidity co- 
efficient of a semi-flexible rod. The procedure will be de- 
scribed in greater detail in the following section. Using 
the effective Hamiltonian H c g t i of Eq. (III.ll), the OSF 
result, Eq. (1.1) is reproduced, with the average density 
of bound counterions as predicted by Manning: 



— 'o + 'e,l 



le,l 



1 



*OSF 



Az 2 k 2 Ib 



(III. 13) 
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IV. BEYOND MEAN FIELD 

The results obtained in the previous section are strictly 
on a mean-field level 0, [45| . In order to go beyond this 
approximation, higher than linear powers of tp h ave to 
be included in the partition func tion, Eq. (ILl). The 
exact partition function, Eq. (tl.l), can be rewritten as 



AH = -an J d 3 r0 m (r) 
(^(r))" + .. 



2! 



(V> b (r))' 



v 
n! 



(IV.l) 



Performing a cumulant expansion 

Z = Tr^^^^e-^ 1 exp [- (AH) c>1 

-Aah 2 ) ^--(ah*) , + . 

2| \ / c,l q| \ /c,l 



3! 



(IV.2) 



where (O n ) c 1 denotes the n-th cumulant. The subscript 
1 indicates that the cumulants are calculated using the 
first-order expansion of AH. For instance, (C 2 ) c j = 

((C 2 ) 1 - (OfX where the moments (O n ) 1 are defined 



according to Eq. ( |III.4| ). 

T he effe ctive interaction neglected by _ff c ff,i of 
Eq. (III. 11) is therefore 



(AH) C<1 I (Ai/ 2 ) c A + i <Atf 3 ) c , + . . . (IV.3) 



3! 



A. Second Order Corrections 



The only term in Eq. (IV.3) which is of second order in 
ijj h is the first term of the first cumulant: 

H 2 = l -an J d 3 r0 m (r) ([^ b (r)] 2 ) i . (IV.4) 

Applying a Fourier transform, Hi can be expressed in 
Fourier space as well 



Hi 



1 f d 3 ki d 3 k 2 , b b , 



(2tt) 3 (2tt) 3 



(IV.5) 



Subs tituti ng ^jj 1 and ipu m t° the correlation expression, 
Eq. ( III. 5 ), yields a correction to the e ffectiv e interaction 
Hamiltonian, 7J e ff,i, obtained in Eq. (pXLll ) 



H 



1 2 n Z 2 



at, 2 



(1 — O2;no) z 



2 a : 

e — k|R(si)— R(s )| e — «|R(s 2 )— H,(s )| 



L L L 

ds y^si /" dS2 

(i 



|R(si)-R(a )| |R(s 2 )-R(so) 
+ z 2 l B n N. 



(IV.6) 



The above result is obtained by using a lower cut-off at 
distance a on one of the three integrations, and expand- 
ing the integral in powers of na, neglecting all but the 
leading term. In the expansion method, tp h represents 
an auxiliary field that is generated by the bound coun- 
terio ns |f27[ Examining the mean- field interaction, 

Eq. (IV.5), each mode of the monomer concentration 
field <fi m interacts with the average of two ip h auxiliary 
fields. The result is a counterion-monomer-countcrion 
interaction, and appropriately, the first term of -ff G ff,2 in 
Eq. (IV.6) has the form of a three-body interaction. As 
explained in the previous section, the chemical potential 
/i does not depend on the polymer conformation R[s). 
As a consequence, the second term of H c g :2 in Eq. ( IV.6 ) 
does not depend on the conformation R(s) and does not 
contribute to the persistence length. 



B. Density of Bound Ions 

Taking into account second-order corrections, the av- 
erage number of bound counterions changes as well 



[Rl 



l 

L 

n 



d 3 (0 b (r)) 5 



n z(l - azn ) — x[R] , (IV. 7) 

a 

where (0) 2 denoted the average of O c alcul ated with the 
second-order Hamiltonian, H 2 of Eq. (IV.5). 

According to the theory of Manning condensation, 
the density of bound counterions for a straight rod is 
n| ot [rod] = n M @ ||. The chemical potential fjt, de- 



fined through no, should therefore satisfy 
ariQ + hqz(1 — azno)^BX[ r °d] = ar 



M 



(IV.8) 



where x[ r °d] is the value of x[R] for the straight rod con- 
formation, in which |R(s)— R(s')| = \s— s' L Su bstituting 
in Eq. ( pV.g| ) the expression for n M , Eq. ([II. 7), yields 



qx[rod] + 1 ± y/(q X [rod] - l) 2 + 4x[rod] 

n ° = — r~T\ .(1V.9J 

In the limit of extremely high salt concentrations, na 3> 
1, all correlation and fluctuation effects are screened out 
completely and decay expo nenti ally with na. The smaller 
of the two solutions of Eq. ( IV.S| ) is therefore chosen in or- 



der that the Manning value nlf* [R] = n l& reproduced 
no[K,a — ^ oo] = n M . As explained in the previous section, 
the chemical potential n (and consequently no = e M /L 
is not expected to depend on the polymer conforma- 
tion. The density of the bound counterions nlf'fR] does, 
howe ver, depend on the polymer conformation through 
Eq. (|IV.7D . 



C. The Persistence Length 



Using the effective Ha milto nian, Eq. ( IV.6 ), and the 
chemical potential, Eq. (IV. 9), it is now possible to re- 



7 



peat Odijk's method and calculate the electrostatic per- 
sistence length. The difference between the effective in- 
teraction energy at a general conformation R(s) as com- 
pared to the straight rod one is 



Heff,i[ rod ] - #69,2 [rod]. 



(IV.10) 



Odijk's method for calculating the persistence length 
requires expanding AH e {{ in small deformations of the 
chain around the straight rod conformation ]30|j . In this 
limit, the distance between points s and s' can be ap- 
proximated as 



|R(s)-R(s')l =s \s-s'\[l-a(s,s')}+0(a 2 ) 



a(s,s') 



1 fs-s' 



24 V P 



(IV. 11) 



where p 3> L is the small overall radius of curvature of the 
fluctuating chain (not to be confused with a spontaneous 
radius of curvature) . The persistence length is then given 
by 

l ca = 2^-AH eS . (IV.12) 
It is, therefore, additive in terms of AH c h : 



l ca = 2^(ff effil [R]- J ff eS , 1 [rod]) + 

+2^ (tf eff , 2 [R] - ff cff , 2 [rod]) . (IV. 13) 

On the other hand, we note that the electrostatic per- 
sistence length is not additive in orders of i[> h (i e ,i 

2 

2^-(if e ff,i[R] — -ffcff,i[rod]) ) as the expression for the 
chemical potential also changes, as compared to the 
mean- field appro xima tion. Inserting in Eq. (IV. 13) the 
value for no, Eq. ( IV.9| ), and expanding the result in pow- 
ers of 1/p the different integrations can be evaluated. 
This requires cutting of all ultra-violet divergences at a 
distance a. In the limit of a < /t' 1 C L we find 



l p — Iq + I, 



c.2 



'OSF 



q(2 



q In na 



(IV. 14) 



where we have expanded l e> 2 in (lnKa)^ 1 and kept the 
two leading terms. Note that we have already taken 
the fi rst or der (mean fie ld) int eraction into account in 
Eq. ( |rV.10D . Equation < \lV.lj ) is our main prediction 
and is depicted in Fig. 1 for different counterion valencies 
z = 1, 2, 3 as a function of na. At low salt concentrations 
(na <C 1) or high q, the persistence length maintains the 
OSF k~ 2 dependence, l c ~ Iosf ~ kT 2 . We find that the 
electrostatic persistence length depends strongly on the 
valency of the counterions. For monovalent counterions, 
l e is usually positive, indicating an effective repulsion be- 
tween the monomers. However, its value is smaller than 



the one predicted by OSF. Introduction of multivalent 
counterions reduces significantly the rigidity of the PE 
and usually l c < 0, indicating an effective attraction be- 
tween monomers. 

The vanishing of the total persistence length l p under 
certain conditions is alluding to the phenomena of PE 
collapse. A full consideration of the rod-globule tran- 
sition requires a more consistent elastic model for the 
polymer chain than the persistence length prescription 
used here. Furthermore, Odijk's method for calculating 
l e does not hold for flexible polymers |2^]. However, the 
condition that total persistence length vanishes, l p = 0, 
is the validity limit of the rod-like regime and is indica- 
tive of some mechanical instability. For instance, using 
parameters applicable to DNA chains: Iq = 500A, 3:1 
salt, a = 1.7A, and Ib = 7 A we get a DNA collapse at 
kT 1 ~ 30A, corresponding to a 3:1 salt concentration of 
about 17m M. 

Expanding the density of the bound ions, Eq. ( [V.Tj ), 
for small deformations of the chain yields 



.M 



1 



In 2 In 



(IV. 15) 



M 



For a straight chain, p — > oo and the Manning value n 
is obtained. However, more counterions are condensed 
on a bent polymer with a finite radius of curvature. This 
enhanced condensation drives a further reduction of the 
persistence length. For a rod- like polymer, p 2 = Ll p /3 
and the correction to n M is of order l/N. This difference 
does, however, have a significant effect on the electro- 
static persistence length, because H c ff,2 is a triple inte- 
gral over the monomers. In order to examine the effect of 
increased condensation, we look at the asymptotic form 
of Eq. ( 1V.14D for q = 1 + Aq, Aq <C 1, in two cases. In 
the first we allow the density of the bound counterions to 
be adjusted according to the equilibri um co ndition with 
the bulk (this is an expansion of Eq. ( |IV.14 ) in power of 
Aq). This is consistent with the general considerations 
of our paper. In the second case we add a constraint that 
fixes the density to be according to the Manning theory 
for all conformations of the polymer. This more restric- 
tive constraint is added in order to make comparisons 
with other models. Expanding in Aq we recalculate l e 
for both cases 

?c = /osf [1 + 0(Ag 2 )] (IV.16) 
C xcd = /osf [1 - [1/ In(/sa)]Ag + 0(Aq 2 )] . 

The lack of a linear term in Aq in the first expression of 
Eq. (IV. 17) indicates that corrections to Manning con- 
densation for bent polymer chains have a substantial 
influence on the persistence length. Fig. 2 depicts the 
excess number of bound counterions as a function of q 
for three salt concentrations (corresponding to na =0.02, 
0.04 and 0.08). 

An interesting effect is charge inversion, where the to- 
tal charge of the polymer with the bound counter-ions 
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changes sign (zn tot < 1/a). The persistence length at 
which charge inversion occurs is given by 



imv 

P ~ a K 2(„ vr ,M 



1 



8n 2 {azn M - 1) In na N 



(IV. 17) 



According to our model, charge inversion will not occur 
on a long (N ^> 1), rod-like polymer, because ip nv oc 1/N. 



D. Higher Order Corrections 

In order to examine the convergence of the loop expan- 
sion used above, we have calculated the next two orders 
of approximation: (ip h ) 3 and (V' b ) 4 - The effective in- 
teractions are obtained by taking into account terms of 
higher orders of if} h m Eq. ( [V.3| ). The 3rd and 4th order 
terms are given as: 



H 3 = 
Hi = 



-an Q I d 3 r0 m (r)^([^ b (r)] 3 ) ^ 



an / d^rfW^^fr))^ 



-50W 



d 3 r^ m (r)L (V- b (r)) 2 



n 2\ 



d 3 r0 m (r)l-(^ b (r))^ 



(IV.18) 



Higher orders of ip h can be taken into account follow- 
ing the same prescription used above for calculating the 
chemical potential and the persistence length for the 
second-order correction. The expression for 4,3 and Z e ,4 
are not detailed here because they are quite complex. 
However, they do not change the polyme r beha vior in 
any qualitative fashion. We find that Eq. ( IV. 14 ), valid 
to second order, accounts for most of the deviation from 
the OSF result. Third and fourth order terms represent 
only a relatively small correction to the second-order one. 
For instance, for q > 10 and na < 0.01, the third-order 
correction is less that 4% of the second-order one. The 
fourth-order correction is again less than 4% than the 
third-order one. 

Furthermore, we find that the convergence is better 
for large q and small na f45}| . This is the more interest- 
ing q regime since, for instance, in DNA solutions with 
trivalent counter-ions we get roughly q ~ 12. 



V. COMPARISON WITH OTHER MODELS 

Our model is closely related to several previous ones. 
Ha and Thirumalai have used a similar loop expansion 
method for calculating the persistence length of polyam- 
pholytes p(| and of bundles of polymers (25) . For the 



case of a single PE, Golestanian et al |35|] took into ac- 
count thermal fluctuations of the bound counterions den- 
sity. This is strictly an all-fluctuations model, which 
is expected to become accurate at high temperatures 
(q = 1 + Aq,Aq <C 1). Their expression for the elec- 
trostatic persistence length is |35fl 



fluct 



h 



Aq 2 n 2 a 2 [l-(q- 1) ln(Ka)]" 



(V.l) 



where we have explicitly omitted terms with a higher 
order dependence on na. 

A second model, suggested by Nguyen et al |j4j as- 
sumes that condensed ions are arranged in a Wigner crys- 
tal, or a strongly correlated liquid on an infinite cylinder 
with diameter d. This picture becomes accurate at low 
temperatures (q ^S> 1). For the case of no-added salt, 
their expression for the electrostatic persistence length 
is: 



I' 



l B V^{d/a) 3/2 . 



(V.2) 



As Nguyen et al did not take into account the effect of 
salt (they considered only counterions), we have slightly 
modified their derivation for the purpose of comparison 
with our model. Equation (13) of Ref. fS] estimates the 
interaction energy of an ion with its Wigner-Seitz cell of 
background charges. Up to a constant of order unity, it 
is found to be 



e(n) 



n 1 / 2 z 2 e 2 



(V.3) 



where n is the average surface charge density on the cylin- 
der. The only change we include in the above equation 
is to assume that in the presence of salt, the interaction 
energy should be proportional to the number of Wigner- 
Seitz cells that reside within a circle with a radius equal, 
roughly, to the screening length kT 1 . Up to a constant 
prcfactor, this modification gives 



n 1 / 2 z 2 e 2 n 



.2 ' 



(V.4) 



Repeating the derivations of Ref. mm with this modifi- 
cation of their Eq. (13), and assuming that d ~ a yields 



OSF" 



(V.5) 



Our loop expansion method can be shown to account 
qualitatively for both limits of the parameter q. As dis- 
cussed in the previous section, we find q to be the rele- 
vant, temperature dependent parameter that determines 
the system behavior. This is the reason we expand the 
results in the two limits g > 1 and q > 1 rather than 
low or high temperatures (k — > 00 or 0). The limits 
of very low or high temperatures are beyond the va- 
lidity range of our model which explicitly assumes that 
L _1 </(< a -1 . However, experimental systems usu- 
ally have a large value of q with a finite screening length. 
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For instance, typical parameter values in experiments 
with DNA segments at room temperature are Ib = 7A, 
a = 1.7 A, K' 1 = 1Q-100A and N > 150. With tri-valent 
counterions (z = 3) we get q ~ 12 ^> 1 and kL > 1. This 
is an example where we can consider the relatively high 
q limit at room temperature. 

In order to compare the three models, we expand the 
expressions for the electrostatic persistence length: ours 

jloop. ;fluct of Ref Q. and jcorr of Ref jjj^ j n thege twQ 

limits. 

Including terms up to linear order in Aq close to q = 1 
we get 



Our expression for the electrostatic persistence, 



?e uct = 'osf [1 - 2 ln( K a)Aq + 0(Aq 2 )] 

C rr =s Iosf + 0(Ag 2 )] 

^ oop = ^OSF [1 + 0(Aq 2 )] . (V.6) 



However, the calculations of Golestanian et al assumes 
that the amount of condensed counterions is according 
to Manning for all polymer conformations. In the loop 
calculation, we have relaxed this assumption and took 
the Manning counterion value only for the completely 
straight rod case. For the sake of comparison, we im- 
pose now this restriction. This has been done already in 
Eq. 0Tfh 



;loop,fixcd = ; ogF ^ _ Aq/HKa) + (Aq 2 )] .(V.7) 

With the new restriction, the linear term in Aq reappears 
but with a different coefficient than in ig° rr of Eq. ( |V.6| ) . 
The differences in the coefficients, as well as the different 
dependence on the cut-off distance a is due to the differ- 
ent methods and approximations used in calculating the 
persistence length. 

The second case is that of large q, for which the models 
give: 



;fluct 

L c 

7Corr 
'e 

^loop 



Iosf [l/[glnM] 2 + 0(l/ g 3 )] 

hsF [~q 2 /V^+0(q)] 



= I 



OSF 



-q 



0(g)] 



(V.8) 



We note that the electrostatic persistence length of the 
two previous models [55| depends very differently on 
each of the fundamental parameters of the system: the 
charge density of the polymer Ib/cl, and the valency of 
the counterions z. Furthermore, their expressions do not 
have similar limits in the two q limits discussed above. 
However, in the limit q > 1, our result is similar to the 
one obtained by the fluctuation model pq |. Conversely, 
in the limit g> 1, our result resembles the one obtained 
by the correlation model |34j . Some discrepancies are ap- 
parent. As explained above, the difference between our 
model and the fluctuations one is mainly due to the dif- 
ferent method used for obtaining the persistence length. 
The difference with the correlations governed model is 
mostly due to the discreteness of the charges assumed in 
Ref. |34f and their specific arrangement in a 2D Wigner 
lattice. 



Eq. (IV. 14), neither vanishes no diverges in the limits 
or low or high temperatures. Fluctuation contributions 
to the electrostatic persistence length vanish in the limit 
T — > 0. Conversely, correlation contributio ns van ish in 
the limit T -> oo. This illustrates that Eq. flrV.14j ) con- 
tains contributions from both fluctuations and correla- 
tions. 



VI. COMPARISON WITH EXPERIMENTS 

Comparison between our expression for the persistence 
length, Eq. ( IV. 14 ), with that obtained in experiments 
0, |[ ||, [l^, |II| is difficult. Although l c .i correctly pre- 
dicts that the persistence length should be smaller than 
OSF, it seems that the reduction we obtain is too large 
as compared with experiments. Actually, at present, we 
are not aware of any other theoretical modeling which 
explains quantitatively the experimental data. 

Measurements of rigid PEs usually involve short DNA 
segments. Experiments show that adding of very small 
amounts of multivalent counterions greatly reduce the 
persistence length of DNA, l p — Iq + l c even bellow its 



bare value Iq |Q, [To| |5(J, indicating a negative electro- 
static persistence length. However, substituting DNA 
parameters (a ~ 1.7A, Iq ~ 500A) and salt concentra- 
tions common to experiments (10A < k^ 1 < IQ OA ) into 



the total persistence length, l p , of Eq. (IV. 14), we find 



that DNA should collapse (l p < 0), for z = 3 or 4. This 
does not agree with experiments where the DNA is still 
in the rigid-rod limit for the same system parameters. 

This discrepancy may be caused by several important 
experimental features which are neglected in our model as 
well as in Refs. |34|, ^5|. DNA segments are prepared in a 
buffer which stabilizes the solution pH and removes free 
divalent calcium ions |30|. The buffer itself contributes 
a finite, non-negligible concentration of monovalent ions. 
In a second stage, multivalent ions such as spermidine 
(z=3) or spermine (z=i) are added. In experiments, the 
monovalent salt concentration may be much higher than 
the multivalent one, making the contribution of the mul- 
tivalent salt to the screening length k" 1 quite negligible. 
In the model we did not take into account mixtures of 
mono- and multi-valent ions. 

According to the Manning-Oosawa model which was 
employed by us, entropy considerations dictate that con- 
densation of multivalent counterions is much favorable 
than monovalent ones. At low concentrations of multiva- 
lent salts, experiments clearly show that this is not always 
the case [|o[ [H], |5^] . At low but finite polymer and mul- 
tivalent salt concentrations (as is usually the case in ex- 
periments) , entropy and the finite size of the counterions 
(spermidine and spermine are relatively large molecules) 
prohibit multivalent counterions from condensing on the 
chain. Generalization of the Manning-Oosawa model to 
account for this effect is not straightforward. 

In experiments the change from a rigid rod-like behav- 
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ior to a flexible one is highly sensitive to the multivalent 
concentration @, 0. The above discussion emphasiz- 
ing the deviation from the Manning-Oosawa model may 
also explain this change-over. At low concentrations of 
multivalent ions, less counterions are condensed than ac- 
cording to the Manning value. Hence, the chain is still 
rigid in disagreement with our prediction. At higher mul- 
tivalent salt and polymer concentrations, where the Man- 
ning value for condensation on a cylinder n M holds, the 
DNA is completely collapsed, making both theory and 
experimental measurements of persistence length useless. 
Therefore, the main difficulty in our (and similar) mod- 
els is the small window of parameter values where the 
model can be applied. Experiments with DNA do not, 
in general, fall in this window due to the strong charging 
of the chains. 

We briefly mention other features not considered in our 
model, and which may influence the persistence length. 
They include the finite size of the counterions [p3[ , the 
ordering of the charges along the polymer chain [ |54| , |55| , 
the concentration profile of the condensed ions around 
the polymer E3| and other, more specific, details of the 
polymer type and ions used in experiments |56|. 



VII. SUMMARY 

We have found significant corrections to the persistence 
length of a single, stiff, strongly charged and long PE, 
as compared to the standard mean-field result of Odijk- 
Skolnick-Fixman. Our method takes into account both 
correlations between condensed ions and thermal fluctu- 
ations. At low salt concentrations, the calculated elec- 
trostatic persistence length Z .2 is proportional to Zosf- 
However, the prefactor, which depends on q = zIb/ci, 
drastically changes the system behavior. For q < 1, 
lc,2 = Iosf is obtained exactly. For 1 < q < 2, the 
electrostatic persistence length l C: 2 is positive, indicating 
an effective repulsion between the monomers. For q > 2, 
the interaction becomes attractive causing a reduction in 
the chain stiffness, Z c ,2 < 0. This observation is in agree- 
ment with the reduction in persistent length observed in 
experiments with multivalent counterions and strongly 
charged polymers. 

We compared o ur res ult for the electrostatic persis- 
tence length, Eq. ( IV. 14 ), with two previous models and 
found that our model takes into account both thermal 
fluctuations and correlations between bound counterions. 
Our model qualitatively agrees with both previous ones 
at different limits of the parameter q. 

It is interesting to note that q = 2 corresponds to the 
case where the average electrostatic interaction between 
bound ions equals k&T. This means that for 1 < q < 2, 
thermal fluctuations are expected to dominate over cor- 



relations. On the other hand, for q > 2, correlations be- 
come more significant. Our conclusion is that although 
thermal fluctuations reduce the (mean field) electrostatic 
repulsion between monomers, they are not sufficient to 
induce effective attraction. In order to correctly describe 
the attractive, collapsed case, correlations between coun- 
terions have to be included. 

We have also obtained the average density of bound 
ions and found that more counterions condense on the 
chain than is predicted by the Manning-Oosawa model. 
The increased condensation has a significant effect on 
the persistence length and cannot be neglected. Fur- 
thermore, we have estimated the conditions under which 
collapse of a rigid PE may occur. The results are reason- 
able and relate, at least qualitatively, to the phenomena 
of DNA condensation. As explained in previous sections, 
our theory cannot yet be directly compared with experi- 
mental results (in particular with DNA). 

We believe that additional work is needed to shed more 
light on the mechanical instability of the chain, indicat- 
ing a rod-globule transition. Different, more complex, 
methods are required in order to calculate the persistence 
length of flexible chains. For instance, the validity range 
of some variational methods are known to be wider than 
that of the OSF theory |27|, |2^], and may better apply 
to flexible chains. Moreover, the entire notion of per- 
sistence length for describing the chain elastic properties 
breaks down close to the instability. The vanishing of the 
persistence length indicates that the method is no longer 
consistent, and higher powers of the radius of local cur- 
vature d 2 R/ds 2 have to be taken into account. 

DNA experiments are usually performed for polymer 
and salt concentrations requiring a more detailed exam- 
ination of the counterion condensation phenomena than 
the simplified Manning-Oosawa model. Some effects that 
are unique to mixtures of different types of ions need to 
be taken into account. A quantitative analysis of this 
phenomenon requires further and rather detailed consid- 
erations and will be presented elsewhere. 
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Figure Captions 



Fig. 1: The electrostatic persistence length l e as func- 
tion of no, according to OSF (dashed line) and our 
of Eq. (IV. 14) (solid line). Valencies are specified next 
to each curve. The parameters chosen are: a = 4A, 
Ib = 7 A, so that q — 1.75z. The negative l c values for 
z = 2,3 indicate a possible collapse transition of the PE 
chain. 



Fig. 2: The excess number of bo und counterions 
L(r4 ot [/»] -n M ) according to Eq. jTV.lBj ), as a function of 
q. The number of excess bound ions is plotted for three 
salt concentrations corresponding to na =0.02, 0.04 and 
0.08. The r adius of curvature was calculated according 
to 4,2, Eq. ([TV.14D, with z = 3, a = 4A and h = 500A. 
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